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Abstract. We investigate the field-angle-dependent zero-energy density of states for YNi2B2C 
with using realistic Fermi surfaces obtained by band calculations. Both the 17th and 18th bands 
are taken into account. For calculating the oscillating density of states, we adopt the Kramer- 
Pesch approximation, which is found to improve accuracy in the oscillation amplitude. We show 
that superconducting gap structure determined by analyzing STM experiments is consistent with 
thermal transport and heat capacity measurements. 



The discovery of the nonmagnetic borocarbide superconductor YN12B2C [I] has attracted 
considerable attention because of the growing evidence for highly anisotropic superconducting 
gap and high superconducting transition temperature 15. 5K. In recent years, Maki et al. [2] 
theoretically suggested that the gap symmetry of this material is s+g wave and the gap 
function has zero points (point nodes) in the momentum space. Motivated by this prediction, 
field-angle dependence of the heat capacity [3] and the thermal conductivity [I] have been 
measured on YN12B2C. The gap symmetry can be deduced from their oscillating behavior. 
Those experimental results 13] were considered to be consistent with the s+g-wave gap. 
However, the present authors [5] recently found that the local density of states (LDOS) around 
a vortex calculated for the s+g- wave gap on an isotropic Fermi surface (FS) is not consistent with 
measurements by scanning tunneling microscopy and spectroscopy (STM/STS) [6j. Therefore, 
we calculated the LDOS around a vortex with the use of a realistic FS of the 17th band obtained 
by a band calculation [7J. We also investigated the density of states (DOS) under zero field 
and the field-angle dependence of the zero-energy DOS (ZEDOS). Consequently, we proposed 
alternative gap structure for YM2B2C and succeeded in reproducing those experimental 
observations consistently [7J. In this paper, we investigate the field-angle-dependent ZEDOS 
taking the 18th band into account in addition to the 17th band considered previously. While 
the so-called Doppler-shift (DS) method was previously utilized in Ref. [7], we adopt here a more 
reliable method [8] for calculating the ZEDOS. 

The DOS is the basis for analyzing physical quantities such as the specific heat and the 
thermal conductivity [9]. For example, the specific heat C/T is proportional to the ZEDOS in 



the zero temperature limit T — > 0. The DOS is obtained from the regular Green function g 
within the quasiclassical theory of superconductivity, which is represented by a parametrization 
with a and b as g = —(1 — ab)/(l + ab). They follow the Riccati equations (h = 1) |10j : 

vp ■ Vci + 2uj n a + aA*a - A = 0, (1) 
vp ■ V6 — 2a) n 6 — bAb + A* = 0. (2) 

Here, t?F is the Fermi velocity, and iuj n = iiv n + (e/c)vp ■ A with the Matsubara frequency oo n 
and the vector potential A. 

To analyze the field-angle-dependent experiments quantitatively, we have developed a method 
on the basis of the Kramer-Pesch approximation (KPA) [8]. This method enables us to take 
account of the vortex-core contribution which is neglected in the DS method. By virtue of it, 
one can achieve quantitative accuracy. We consider a single vortex situated at the origin of the 
coordinates. To take account of the contributions of a vortex core, in the KPA we expand the 
Riccati equations (pQ) and ([2]) up to first order in the impact parameter y around a vortex and 
the energy oj n [5] (ioj n — > E + irj and y is the coordinate along the direction perpendicular to 
v-p). By this expansion, one can obtain an analytic solution of the Riccati equations around a 
vortex. We then find the expression for the angular-resolved DOS [8] 

N(E,a M ,0 M ) = ™( [ ^^f} % ) . (3) 
2ir 2 £o\J \v F \ (E - E y ) 2 + rf / gp 

Here, the azimuthal (polar) angle of the magnetic field H is om (#m) hi a spherical coordinate 
frame fixed to crystal axes, and dSp is an area element on the FS [e.g., dSp = k F sax0d(j)d6 for a 
spherical FS in the spherical coordinates (k,4>,9), and dSp = kp a bd<j)dk c for a cylindrical FS in 
the cylindrical coordinates (kab,cj),k c )]. In the cylindrical coordinate frame (r,a,z) with z || H 
in the real space, the pair potential is A = AoA(fcF) tanh(r/£o) exp(za) around a vortex, Ao 
is the maximum pair amplitude in the bulk, A = |A|, and (• • -)sp = Jq 11 rdr J Q 27r • • • da/ (vrr^) is 
the real-space average around a vortex, where r a /£o = \J H C 2 / H [H C 2 = $o/(7r£o)> $o = 7Tr a^]- 
x = rcos(a — 6 V ), y = rsin(a — 6 V ), and E y = AoA 2 y / '(£,oh). 9 v (kp,aM,0M) is the angle 
of Upj_ in the plane of z — 0, whore ot a/iid. V are measured from a common axis (5j |7j. 
vp± is the vector component of vp(kp) projected onto the plane normal to H = (ckmj^m)- 
\vp±(kp,au,du)\ = ^fo(«m, #m)M^f, «m, ^m) and vpo is the FS average of \v F ±\ [7]. is 
defined as £o = ^Fo/( 7r Ao). We consider here a clean SC in the type-II limit. The impurity 
effect can be incorporated through the smearing factor i]. 

We calculate the angular dependence of the ZEDOS (E = in Eq. ([3])) for YM2B2C with 
using the band structure calculated by Yamauchi et al. When integrating Eq. (J3]), the 

band structure is reflected in dSp and vp (vp±, vpo, and h are obtained from vp). In YM2B2C 
there are three bands crossing the Fermi level, which are called the 17th, 18th, and 19th bands. 
The electrons of the 17th band predominantly contribute to the superconductivity, since the 
DOS at the Fermi level for the 17th, 18th, and 19th bands are 48.64, 7.88, and 0.38 states/Ry, 
respectively [11]. We consider here the FSs of the 17th and 18th bands and neglect the 19th 
band because the DOS is substantially small on the 19th FS. We use an anisotropic s-wave 
gap structure obtained in our preceding paper (Eq. (27) in Ref. [7|) for the 17th FS. This gap 
structure was determined by analyzing STM observations [6j. As for the 18th FS, we assume 
an isotropic gap with an amplitude equal to the maximum gap on the 17th FS, according to 
angular-resolved photo emission measurements [12j . 

We rotate the applied magnetic field in the basal plane perpendicular to the c axis (9m = 7r /2) 
and investigate the azimuthal angle «m dependence of the ZEDOS. First, we show the partial 
ZEDOS on the 17th FS only. As seen in Fig. [IJa), the oscillation amplitude is of the order 9 %. 



(a) DOS with the 17th band 



(b) DOS with the 18th band 



(c) DOS with the 1 7th and 1 8th bands 
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Figure 1. Angular dependence of the ZEDOS for YNi2E>2C (a):with use of the 17th band, (b): 
that of 18th band and (c) that of the 17th and 18th bands. The magnetic field tilts from the c 
axis by polar angle 9m = 7r /2- f] = 0.05Aq and r a = 7£o- 



Previously we have calculated the ZEDOS by the DS method under the same condition and 
obtained the amplitude of the order 30 % [7j, which was too large in comparison to experimental 
results 0H]. The KPA adopted in this paper yields substantial improvement in the amplitude. 

Second, we show the partial ZEDOS on the 18th FS only. As seen in Fig. (U[b), the cusp-like 
minima appear. The origin of the cusp-like structure is different from that of the 17th-FS case, 
since the gap is isotropic on the 18th FS while it is anisotropic on the 17th FS. The cusp-like 
minima are due to the FS anisotropy on the 18th FS. According to the band calculation (Fig. 
2(b) in Ref. [IT]), the shape of the 18th FS is like a square as viewed from the k c axis in the 
momentum space. The sides of this square are parallel to the k a or kb axis. Therefore, the region 
where the Fermi velocity is parallel to the k a or kf, axis has a high proportion of the area on the 
18th FS. Now, the quasiparticles with the Fermi velocity parallel to the magnetic field do not 
contribute to the ZEDOS (e.g., see Ref. [8]). As a result, when the magnetic field is rotated near 
the directions parallel to those axes, the change in the ZEDOS is drastic, leading to cusp-like 
minima. Udagawa et al. [13] previously obtained cusp-like minima in the field-angle-dependent 
ZEDOS by using a square-like FS model. The origin of their cusp-like minima is the same as 
that of the present result for the 18th FS. It should be noted that the DOS ratio r of the 18th 
FS to the 17th FS is r = N F18 /N F17 = 7.88/48.64 ~ 0.16 [IT]. The contribution of the 18th FS 
is not dominant, and therefore the cusp-like minima observed in the experiments [Hill] probably 
cannot be attributed to the square shape of the FS only. 

Finally, we show the total ZEDOS on both the 17th and 18th FSs. As seen in Fig. Hfc), the 
oscillation amplitude is of the order 8 % and the overall behavior is almost similar to that for 
the 17th FS shown in Fig. Ufa). In the thermal conductivity measurements in YM2B2C, the 
oscillation amplitude in the field 1 T and 0.5 T is of the order 2 % [1] and 4.5 % [H], respectively, 
at the same temperature 0.43 K. In the heat capacity measurements [3], the oscillation amplitude 
in 1 T at 2 K is of the order 4.7 %. The KPA is an approximation appropriate in the limit of 
the low temperature and low magnetic field. Therefore, in terms of the order of the oscillation 
amplitude, our result could be considered to be consistent with that of those experiments. As 
for the directions of the minima, the result is also consistent with the experiments. Those mean 
that the assumed gap structure determined by analyzing STM experiments (Eq. (27) in Ref. 
[7]) is consistent with the heat capacity [3] and the thermal transport [I] observations. So far 
we have found here that inclusion of the 18th-FS contribution does not drastically alter our 
previous conclusion [TJ, and the use of the KPA improves the value of the amplitude compared 
to our result [7J by the DS method. We should note, however, that the sharpness of the cusp 
is rather weak in Fig. (U[c) in comparison to the experimental observations 0]. This point 
may be resolved by considering carefully the k c dependence of the gap structure. The gap 



structure used here was determined to explain the azimuthal dependence of the observations 
in the basal plane [7J. A modification of the gap structure would be necessary to explain the 
observed polar-angle dependence [3]. This is left for future studies. Here, it should be noted that 
such a detailed discussion is almost impossible till the KPA [8] is developed as a new method 
with high precision. 

In conclusion, we calculated the field-angle-dependent ZEDOS on the FSs from 17th and 
18th band of YM2B2C. We adopted the gap structure [7J determined by comparison with STM 
observations [BJ. The oscillation amplitude was found to be of the order 8 %. This amplitude does 
not deviate so much from the heat capacity [3] and the thermal transport [HE] measurements. 
As for the directions of the minima, the result is also consistent with those experiments. This 
suggests that the assumed gap structure on the 17th FS is essentially appropriate for the in-plane 
properties. Further modification of the gap structure is expected in terms of the k c dependence. 
The use of the KPA improved the value of the amplitude compared to our previous result [7] by 
the DS method. The KPA [8] appears to be an efficient method with high precision for analyzing 
experimental data. 
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